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Abstract 



\ We present a method of treating the interaction of a single three-level ion with 

two laser beams. The idea is to apply a unitary transformation such that the exact 
transformed Hamiltonian has one of the three levels decoupled for all values of the 
detunings. When one takes into account damping, the evolution of the system is 
governed by a master equation usually obtained via adiabatic approximation under 
the assumption of far-detuned lasers. To go around the drawbacks of this technique, 
we use the same unitary transformation to get an effective master equation. 

Ti.: 

^ ' 1 Introduction 

> ' 

In recent years, much progress has been made in cooling and trapping tech- 
niques for neutral atoms and for ions [1,2,3,4,5,6,7]. In fact, in a trapped ion the 
center-of-mass dynamics gets entangled with the internal degrees of freedom 
and some peculiar aspects of the vibrational response have been successfully 
exploited for experimentally generating Fock [8], coherent and squeezed [9], 
and Schrodinger catlike states [10], proposing theoretical schemes for engineer- 
ing several nonclassical states [11,12,13,14,15,16,17,18,19,20], realizing tomo- 
graphic reconstructions of the density matrix [21,22,23,24], and characterizing 
a variety of quantum effects [25,26,27,28,29]. This is of fundamental inter- 
est, since it brings to the forefront issues involving the relationship between 
quantum and classical physics, but also offers potential applications for, e.g., 
precision spectroscopy [30] or quantum computation [31,32,33,34]. 

In modeling typical experiments one considers a three-level atomic system 
interacting with two laser fields (Raman scheme [35,36]) and reduces it to a 
two-level problem on the assumption of large detunings by using the adiabatic 
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elimination [37]: the effective Hamiltonian obtained in this way has the form 
of the usual Jaynes-Cummings model. 

Adiabatic elimination has been criticized on several grounds [38,39,40,41], and 
other methods of deriving effective Hamiltonians exist [42,43]. In this spirit, 
we have recently proposed an alternative approach that involves using a uni- 
tary transformation (in fact, a nonlinear rotation) to obtain an equivalent 
Hamiltonian for which one level decouples [44,45]. The transformation can 
be exactly found and gives the same results as the adiabatic elimination (ex- 
cept for including intensity-dependent Stark shifts) when it is evaluated up to 
second-order terms in coupling constants. 

To take into account the effects of damping in Raman cooling schemes, the 
standard way of proceeding is to start from the master equation for the three- 
level system and adiabatically eliminate the far-off resonant level. The details 
are described in many different text books [46,47]. Unfortunately, it is known 
that this treatment is not valid in many regimes of physical interest and other 
approximations are required [48,49,50,51,52,53]. 

The main purpose of this paper is to show how our approach of nonlinear 
rotations allows one to go around these drawbacks in a natural way. Our 
strategy can be stated in very simple terms: starting from the exact master 
equation for the three-level model, we apply to it the same unitary transfor- 
mation leading to the effective Hamiltonian, obtaining in this way what we 
call an effective master equation [54]. Here we fully investigate this approach 
and present numerical evidences of its validity. 



2 Physical system and model Hamiltonian 

In the interest of retaining as much clarity as possible, we first recall some well- 
known facts [55,56] about the system we wish to treat here, which consists of a 
three-level trapped ion in the A configuration with energy levels E < E 1 < E 2 , 
as shown in figure 1. As usual, to describe this system we use the operators 

Sv = \m, (i) 

where \i) denotes the eigenstate of the iih atomic level. One can easily check 
that they satisfy 

[Sij, Sm] = 5jkSu — SuSkj , (2) 

which correspond to the commutation relations of the algebra u(3). Obviously, 
the three "diagonal" operators Su measure level populations, while the "off- 
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Fig. 1. Energy scheme of a three-level A- type ion interacting with two laser fields, 
coupling the two ground states to a common excited atomic state via a Raman 
transition. Here we illustrate the case of driving the first blue sideband. 



The levels |0) and |1) are metastable and coupled by stimulated Raman tran- 
sitions via two classical optical fields (of frequencies u> a and Uf,) of the form 



where from now on the index £ runs the values a and 6, x is the position 
operator associated with the center-of-mass motion and (ft is the phase of the 
laser field i at the mean position of the ion. 

The ion is trapped in a harmonic potential. Therefore, the center-of-mass 
motion can be described in terms of annihilation and creation operators of 
vibrational quanta (phonons) in the usual way 



where v q represents the oscillatory frequency along the gth direction, M is the 
ion mass, and Ax q is the width of the ground-state wave function. 

The Hamiltonian that describes the system is H = H cm + H{ on + V, where 



diagonal" Si* generate transitions from level i to level j. 




(3) 




(4) 




4m = fr u <i ®l 



(5) 



V = h[g a {±) e-^&a + (£(*) e*""*&o] 
+ % 6 (x) e-^Sn + # 6 *(x) e^'&i] . 
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The interaction term V is written in the rotating-wave approximation and the 
coupling constants are 

g e (x) = K, e £ e exp[i(k £ • x - (p e )] , (6) 

where Kg is the corresponding dipole matrix element in the direction of the 
driving field. 

Because one has the trivial constraint S 00 + Su+ S 22 = 1, only two populations 
can vary independently. Eliminating the population of the level \2) we can 
recast the Hamiltonian as H = H + Hi nt , with 

H = Y^ frv q a\a q - h(A a + uj a )S o - ^(A 6 + u b )Su , 
q 

(7) 

H int = h[g a {±) e-^Sv + g* a (x) e^S 20 ] 
+ % 6 (x) e-^'Su+gltk) e^S 21 ], 

where we have defined the following detunings 

hA a = E 2 -E - hu a , HA b = E 2 -E 1 -hw b . (8) 

This Hamiltonian contains terms oscillating rapidly in time at frequencies u a 
and u b , which can be eliminated by going to a rotating frame. The final result 
is 

H o = Yl hv q a\a q - h(A a S 00 + A b S u ) , 
1 

(9) 

H iQt = h[g a (x)S 02 + g*(x)S 20 } + H[g b (x)S 12 + g* b {x)S 21 ] . 
This is the basic Hamiltonian that will be used in our subsequent analysis. 



3 Effective Raman Hamiltonian in the dispersive limit 

The standard treatment assumes that level |2) is far off resonance and proceeds 
via adiabatic elimination to obtain an effective two-level Raman interaction 
Hamiltonian with an intensity-dependent coupling between levels |0) and |1). 
In Refs. [44] and [45] we have claimed that such a procedure has drawbacks 
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and proposed instead an alternative technique involving nonlinear rotations. 
For the problem at hand, we introduce the unitary transformation 



U = TEf\ (10) 

where T = exp[J(x)] and 

J(x) = [e a (±)S 02 - e* a (±)S 20 ] + [e b (±)S 12 - e* b (±)S 21 ] . (11) 

Here the parameters £ a (x) and £ b (~x.) are defined by 

= (12) 



To interpret this operator T, we note that (S 02 , S 20 ) and (Si 2 , S 2 i) are raising 
and lowering su(2) operators that correspond to the allowed transitions 0^2 
and 1^2. However, these two dipoles are not independent, since Eq. (2) 
imposes highly nontrivial couplings between them. In consequence, T can be 
seen as a "rotation" acting on the subspace of these two dipoles. 

By using the well-known expression 

e^e-^E" (13) 

n=0 n - 



where = [A, B^^ 1 " 1 ], and B^ = B, the exact transformation law (10) 
has been found in Refs. [57] and [58] (see also [59]). The explicit expression 
is complicated although is valid for any values of the detunings. Since we are 
assuming that level |2) is far off resonance and, consequently, the ratios £ a (x) 
and £&(x) can be taken as small quantities, the series (13) can be evaluated 
keeping only up to second-order terms. By applying this to (9) we finally 
obtain an effective Hamiltonian H = Hq + 'Hint, where 



H = £ H a\a q - hU a + S m -h(A b + S u , 



1 



(14) 



7i in t = -- + ^- ) [.y ( ,(x)<y*(x)6' (1 ;i -f .y;;(x)y,,(x),S' in ] . 



In the usual experiments the wave vector difference is chosen to be parallel 
with the X direction of the trap, so (k a — k b ) • x = 8k x and the interaction 
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couples only the motion in X direction to the internal state of the trapped 
ion. Using the effective inversion between levels |1) and |0): 



we can write 



H = hv a) a + hS , (16a) 
Ml 



Hint = - — (e idkx S Q1 + e~ ldkx S 10 ) , (16b) 



with 



In terms of the phonon raising and lowering operators we rewrite (16b) as 

^nt = ^[e">( a+at ) S 01 + e-^ a+at ) S w ] , (18) 



where the Lamb-Dicke parameter is defined as 



"Wass;- (19) 



and represents the ratio between the recoil energy and the quantum vibra- 
tional energy, both taken in the X direction. It is worth observing that the 
second-order corrections to this effective Hamiltonian vanish, so (18) accu- 
rately describes the system dynamics up to times t ^ h/{g^£f). 

In the interaction picture relative to Tl we finally get 

U int = -^{eWa+a*)-**] S 01 + e -^(a + ^)-5t] § w} (2Q) 
2 



By tuning the frequency difference 5 to an integer multiple of the trap fre- 
quency v, 5 — {n' — n)u, we can resonantly drive transitions from |0, n) to 
|1, n'), where the ket \j, n) indicates the nth vibrational Fock state in the elec- 
tronic state j (j = 0, 1). In this case, 7i int is dominated by a single stationary 
term. The exponent exp[ii](a + al)] in Eq. (20) contains all powers of a and 
a*. However, all contributions with m ^ n' — n oscillate rapidly and average 
out when v is much larger than Q. We assume the Lamb-Dicke limit, in which 
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r\\Jn + 1 <C 1. In the relevant case of the first red sideband 5 = v, we get to 
lowest order in r] 

Ml 

'Hmt = -iv-^(aSoi - wSio) , (21) 

which is the familiar Jaynes-Cummings Hamiltonian. Similarly, there is a first 
blue sideband, corresponding to an anti- Jaynes-Cummings Hamiltonian and 
higher-order sidebands, but in the rest of this paper we shall be mainly con- 
cerned with the model Hamiltonian (21). 



4 Damping in terms of an effective master equation 

For many purposes, the coherent control of the vibrational dynamics plays a 
crucial role. For the model discussed in the previous section damping effects 
have been observed that even occurred under almost ideal conditions [60]. 

The presence of the two lasers causes the appearance of a coupling between 
internal (electronic) and external (center of mass) degrees of freedom of the 
trapped ion. However, when the Lamb-Dicke parameter is very small, the dy- 
namics due to this coupling is slow compared with the internal dynamics that 
may be adiabatically eliminated. This reduction leads to a master equation for 
the motional degrees of freedom, where the involved transition rates depend 
on steady-state expectation values of internal operators [37]. 

To take full advantage of the method outlined in the previous section, we start 
from the density operator p for the external and internal degrees of freedom 
of the three-level ion and make the hypothesis that the dynamics is described 
by a general master equation of the Lindblad type [61]: 

j t P = ^[H, P ] + ^C[S 20 ] p+^C[S 21 ] p, (22) 

where 7 a and 7& represent the decoherence rates for the processes associated 
with the coupling of the dipoles with a zero-temperature bath and C[C] is the 
Lindblad superoperator 

C[C] p = 2CpC ] - {C^C, p} . (23) 

Note that (22) describes an irreversible evolution of the system at different 
rates for each channel. It is implicitly assumed that the both dipole moments 
are orthogonal to each other; that is, d 2 • d^ 2 = 0, where djj are the transition 
dipole matrix elements [62]. 
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It seems natural to ask how this equation is transformed by the same unitary 
operator T leading to the effective Hamiltonian in (10). Let us denote the 
effective density matrix by 

g = fpf ] . (24) 
Taking into account that S 20 p = S 2 ip = 0, we get up to second-order terms 
fS Q2 f ] = S 02 - e* a (x)S 00 - e* b (x)S 01 , 

(25) 

TS 12 f j = S 12 - e* b (x)S u - e* a (x)S 10 . 
Then, if we apply T to (22), we obtain the effective master equation 

^-,Q= 4[^int, q] 

at in 

+ y>C[e a (£)S'oo + e b (x)S 10 ] g 

+ ^C[e b (x)S n +e a (x)Soi} g, (26) 

where g is expressed in the interaction picture. 

In the Lamb-Dicke limit we can make the approximation 

e*(x) ~e t (l + ritx), (27) 

with [compare with (19)] 



* = Wdb- (28) 



and the master equation takes the simpler form 

^g= [Hint: Q] 

+ ^^C[S 00 ] g+l^ljC^} g 
+ ^ — C[S 10 ] g + — ^—C[S 01 ] g 

+ ^lC[e a (x)S 00 ,e b (x)S 10 ] g + ^JC[e b (x)S u , e a (x)S 01 ] g. (29) 
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This equation has a very suggestive and transparent physical meaning: the 
terms £[5*00] Q and g describe pure phase dissipation, meanwhile the 

terms £[6*10] Q and £[Soi] Q describe dissipative-like transitions from the level 
|0) to |1) and back, which leads to a stationary distribution of population in 
these levels (incoherent mixing). These two terms seem to simulate the cou- 
pling to a thermal bath at finite temperature, but there is a strong difference: 
now the corresponding rates for the processes from level |0) to |1) and from 
level |1) to |0) are not related by a Maxwell-Boltzmann factor of the form 
nj (n + 1). All these terms have a purely atomic nature and appear in a simi- 
lar master equation derived in a different context by Di Fidio and Vogel [63], 
who interpreted them in terms of quantum jumps. 

Besides, we have also "crossed terms" described by superoperators /C, which 
contribute substantially (i.e., the corresponding terms are time independent 
in the rotating frame) only when some resonance conditions discussed after 
Eq. (20) hold. For the first red sideband we are considering here, they take 
the form 



K,[e a (x)SoQ,e b (x)S 10 \ g = ie a e* b 2(rj a SooagS i - r) b S 00 gaS 01 ) 

- r](aS iSoog + gaS 01 S 00 ) + h. c. , 

(30) 

/C[£ 6 (£)S'ii,£ a (£)S'oi] g = ie b e* a 2(i] b S 11 a^ gS w - i] a S 11 ga ] S w ) 

+ r](a'S 10 S 11 g + ga^S w Sn) + h. c. , 

where we have retained only the dominant terms in the parameter r\ = r\ a — r\ b . 

To test our theory, we have numerically integrated the master equation (29) 
using the Quantum Optics Toolbox produced by S. M. Tan [64]. In the typical 
experiments at NIST [22], a single 9 Be + ion is stored in a RF Paul trap with 
a secular frequency along X of u/2ti ~ 11.2 MHz, providing a spread of the 
ground state wave function of Ax ~ 7 nm, with a Lamb-Dicke parameter of 
rj ~ 0.202. The two laser beams, with 0.5 W in each one, are approximately 
detuned A/2tt ~ 12 GHz, so that Q/2n ~ 475 kHz. With these data we find 
e a ~ e b ~ 0.01, so they can be considered as small parameters, as assumed in 
the previous Section. We take also + 7b)/27r ~ 19.4 MHz (which is about 
one linewidth of the transition) and 7 a = 7^. The observable measured in all 
these experiments is the fluorescence signal, which is the probability -Pj.(t) of 
occupation of the electronic level |0). This probability may be written as 

Pl(t) =£<0,n|^)|0,n). (31) 



Let us consider first the case in which the ion starts in a Fock state with 



9 




10 



20 30 
t (us) 



40 



50 



Fig. 2. Evolution of the P[{t) for an initial Fock state with no = 1 driven by a 
Raman interaction. The parameters are the same as in the experiment [8]: v/2n ~ 
11.2 MHz, A/2vr ~ 12 GHz, 77 ~ 0.202, n/2ir ~ 475 kHz, and (j a + 7&)/2vr ~ 
19.4 MHz with j a = -y b . 

n excitations. The result for n = 1 appears in Fig. 2, showing clearly the 
existence of damped Rabi oscillations. We have numerically checked that in 
this case the role played by the "crossed terms" (30) is insignificant. Similar 
results have been found in Ref. [63] by solving the master equation with quan- 
tum trajectory methods [65,66,67]. The advantage of this approach is that the 
damping can be understood without the need of introducing phenomenolog- 
ical noise. Stochastic models [68,69], leading to a dispersivelike decoherence 
dynamics [70], have also been used for the same reason. We stress that our the- 
ory gives essentially the same results by resorting only to two pure Lindblad 
terms of very easy interpretation. 

An intriguing result found in the experiments of Ref. [8] is that the fluorescence 
signal, for initial Fock states, may be approximately modeled by 



P^t) ~ l - [l + cos(2f2 no t)e-W' 



(32) 



where Vt no is the associated Rabi frequency and the phenomenological decay 
constants 7„ were fitted as 7„ ~ 70(^0 + l) ' 7 . This exponential decay can 
be inferred with good accuracy from a numerical analysis of our simulated 
data [70]. 

In Fig. 3 we show our results for the evolution of the fluorescence signal for 
an initial coherent state with an average vibrational number n — 3. The 
graphic reproduces all the salient features of the experiment [22], although 
for a perfect fitting a more precise value of our parameters r\ a and r/& would 
be needed. The influence of the terms (30) is again very small. 

In Fig. 4 we have plotted the dynamics of the same coherent state with n = 3, 
but in a larger time scale and for two different values of ja/lb- We see that 
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Fig. 3. Evolution of the P[{t) for an initial coherent state with an average vibrational 
number n = 3 driven by a Raman interaction. The rest of the parameters are the 
same as in Fig. 2. 



the stationary limit of the oscillations is not 0.5, due to the phenomenon of 
incoherent mixing mentioned above: while the details of the collapse are almost 
insensitive to the values of the ratio 7 a /7&, as time goes by the quasi-stationary 
values of the population tend to be different. This cannot be reproduced by 
using the standard master equation approach [20]. The "crossed terms" are 
important here in preventing that for larger values of 7a/7& the solutions would 
differ too much. 

To confirm that this is, in fact, an incoherent effect that affects only to the 
population dynamics, in Fig. 5 we show the evolution of the coherence £oi 
for the same initial state, and we clearly see that the values of •y a /lb do not 
influence at all to the dynamics at the large of the dipole moment. 



1.00 




0.25 - 

0.00 -J 1 1 1 1 

1 2 3 4 5 

t (ms) 

Fig. 4. Evolution of the P\(t) for the same initial coherent state with an average 
vibrational number re = 3 as in Fig. 3 in a different time scale. The solid line 
represents the case J a /jb = 1 ; while the dashed line is for j a /^b = 3. 
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Fig. 5. Evolution of the coherences (goi) for the state as in Fig. 4. 
5 Concluding remarks 

What we expect to have accomplished in this paper is to present a comprehen- 
sive method of treating two-photon stimulated Raman transitions in a single 
trapped ion. Our approach is based on the application of a unitary trans- 
formation and leads in a natural way to an effective master equation of the 
Lindblad type with a clear physical interpretation. 

In the framework of this description, we have shown that the evolution of 
populations and dipole moments predicted by our theory is in good agreement 
with the data of realistic experiments. The appearance of an incoherent mixing 
that leads to a stationary redistribution of populations has been also studied. 
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